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We review our recent works about "pasta" structures following the first- 
order phase transition in dense matter, which correspond to the structured 
mixed phases with geometrical symmetries. Three kinds of phase transitions at 
different density ranges are examined as the stages of pasta structures: liquid- 
gas phase transition at subnuclcar density, kaon condensation and hadron-quark 
phase transition at high density. Charge density as well as particle density is 
non-uniform there. A consistent treatment of the Coulomb potential and the 
particle densities is presented and a peculiar role of the Coulomb potential is 
elucidated: the physical picture of the Maxwell construction will be effectively 
recovered. It largely influences the density regime of pasta structures by the 
charge screening effect. 
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1 Introduction 



Recently many efforts have been made to reveal "new form of matter" under 
extreme conditions such as high-temperature and/or high-density. We may find 
such form through high-energy heavy-ion experiments or observation of compact 
stars. In this review article we discuss non-uniform structures of cold matter at 
high density, which is relevant for compact stars. Actually compact stars are 
good laboratories for nuclear physics and elementary particle physics to study 
matter at extreme conditions |2 El EH • 

Nuclear matter in the ground state, which is approximately realized in 
atomic nuclei, consists of the same number of protons and neutrons and has 
the density of po ~ 0.16 fm~ 3 "the normal nuclear density" or "the saturation 
density" , and the binding energy of around —16 MeV. This is called "saturation 
property" of nuclear matter. However, matter in the stellar objects has variety 
of density and chemical component due to the presence of gravity. In fact, at 
the surface of neutron stars, there exists a region where the density is lower than 
po over a couple of hundreds meters. The pressure of such matter is retained by 
degenerate electrons, while baryons are clusterized and have little contribution 
to the pressure. Due to the gravity the pressure and the density increase in 
the inner region (in fact, the density at the center amounts to several times 
the normal nuclear density). Charge neutral matter consists of neutrons and 
the same number of protons and electrons under chemical equilibrium. Since 
the kinetic energy of degenerate electrons is much higher than that of baryons, 
the electron fraction (or the proton one) decreases with increase of density and 
thus neutrons become the main component and seep out from the clusters. In 
this way baryons come to contribute to the pressure as well as electrons. At 
a certain density, other components such as hyperons and strange mesons may 
emerge. For example, negative kaon condensation, expected to be of a first-order 
phase transition (FOPT), remarkably softens the equation of state of matter. 
At even higher density, hadron-quark deconfinement transition may occur and 
quarks in hadrons are liberated. This phase transition is also considered to be 
of first-order. 

If we assume uniform matter, a FOPT brings about a thermodynamical 
instability. In other words, matter should have the mixed phase around the 
critical density. Consider isotherm of water at room temperature for example. 
At low density (or pressure) water exists in the vapor phase. With increase of 
density (or pressure), droplets of water appear in the vapor phase. Finally all 
the volume is filled with liquid water at a certain density (or pressure). The 
equation of state (EOS), i.e. the relation of pressure and density, at low density 
and high density is given by those of the vapor and liquid phases, respectively. 
In between, the EOS of the mixed phase is obtained by the Maxwell construc- 
tion. Required conditions are equalities of chemical potential, temperature and 
pressure between the vapor and liquid phases. 

Phase transitions in nuclear matter are different from that in water. First, 
nuclear matter consists of two chemically independent components, i.e. baryons 
and electrons. Then the equalities of both baryon and electron chemical po- 
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Figure 1: Schematic picture of pasta structures. Phase transition from blue 
phase (left-bottom) to red phase (right-bottom) is considered. 

tcntials between two phases are required by the Gibbs conditions in the mixed 
phase. Therefore the EOS of the mixed phase cannot be obtained simply by the 
Maxwell construction, which is relevant only for single component. Secondly 
those components are electrically charged and the mixed phase is no more uni- 
form. This point is important to the geometrical structure of the mixed phase. 
In the case of the mixed phase of electrically neutral particles, each phase has 
arbitrary geometry when surface tension is negligible or two phases are sepa- 
rated when the surface tension is strong. On the other hand the mixed phase in 
nuclear matter should possess some specific size and shape due to the balance 
between the surface tension and the Coulomb interaction. To minimize the sur- 
face energy plus the Coulomb energy matter is expected to form a structured 
mixed phase, i.e. a lattice of lumps of a phase with a geometrical symmetry 
embedded in the other phase. 

At very low densities, nuclei in matter are expected to form the Coulomb 
lattice embedded in the electron sea, that minimizes the Coulomb interaction 
energy. With increase of density, "nuclear pasta" structures (see Fig. ^) emerge 
as a structured mixed phase [5] in the liquid-gas phase transition, where stable 
nuclear shape may change from droplet to rod, slab, tube, and to bubble. Pasta 
nuclei are eventually dissolved into uniform matter at a certain nucleon density 
below the saturation density po — 0.16 fm~ 3 . The name "pasta" comes from 
rod and slab structures figuratively spoken as "spaghetti" and "lasagna" . Such 
low-density nuclear matter exists in the collapsing stage of supernovae and in 
the crust of neutron stars. Supernova matter is relevant to liquid-gas transition 
of non-beta-equilibrium nuclear matter with a fixed proton mixing ratio and 
the low-density neutron star matter is relevant to neutron-drip transition of 
beta-equilibrium nuclear matter. 

The appearance of the pasta structures could have important consequences 
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for the various neutron star phenomena or supernova. If they appear at an early 
stage of the neutron star evolution nuclear pasta structures could cause a drastic 
change of the neutrino opacity and consequently influence the subsequent 
neutron star cooling [71 [HIE]- It may also affect the matter resistance to the 
stress and consequently the glitch phenomena JUL etc. 

A number of authors have investigated the low-density nuclear matter using 
various models [3 HH E2J H3I HU [TSJ [TBI HZl EBI EU- Roughly speaking, the 
favorable nuclear shape is determined by a balance between the surface and 
Coulomb energies. In most of the previous studies the rearrangement effect of 
the density profile of the charged particles due to the Coulomb interaction has 
been discarded. In Ref. EU the electron screening effect has been studied and it 
has been found that this effect is of minor importance. However, the rearrange- 
ment of the proton profiles as the consequence of the Coulomb repulsion was not 
shown up in their model. The structured mixed phase is also expected in the 
phase transitions at higher density, like kaon condensation and hadron-quark 
phase transition. In these cases, the charge screening effect may be important 
because the local charge density can be significantly high. 

As it was argued in a number of works, the FOPT to a K~ condensate state 
might occur in neutron stars at densities several times larger than the nuclear 
saturation density p (see EH 1221 E31 12H HSI EHl 03 HH1 123 EDI EH E21 E31 
and refs therein). Refs E2I ESI EH E3 I2H ESI EOI EH E2] and many 
other studied possibility of the s-wave kaon condensation, whereas refs |251 E31 
studied different possibilities for the p-wave condensation. It was concluded 
that the kaon condensation in neutron stars may indeed occur and it is the 
FOPT. Such phenomena would lead to interesting consequences in the physics 
of neutron stars: softening of EOS may give a possibility of the delayed collapse 
of protoneutron stars to the low-mass black holes, and the nucleon Urea process 
under background kaons may give a fast cooling mechanism of neutron stars 

mEHmE3E§IE3E3ES]. Ka0n condensation as a FOPT may also lead 
to an expected existence of a wide region of the mixed phase, as was suggested 
in earlier works, cf. |281 1291 15U| . which consists of kaon condensed phase (dense 
phase) and usual nuclear matter (dilute phase) with some geometrical structures 
(kaonic pasta). 

Next we consider the hadron-quark mixed phase. Within the currently ac- 
cepted theory, the deconfincmcnt transition is believed to occur in hot and/or 
high-density matter, although its mechanism is not yet well understood. Many 
people have studied this transition with model calculations and first-principle 
calculations, like lattice QCD 0U|. Above the critical temperature and/or den- 
sity, the quarks inside each hadron are liberated, and we can understand the 
resulting situation as quark matter interacting through gluon exchange. Theo- 
retically, static and dynamic properties of quark matter have been extensively 
studied for quark-gluon plasmas (QGP), color superconductivity ^[021 an d 
magnetism 031 SI I45| . Experimentally, evidence of quark matter has been 
searched for in relativistic heavy-ion collisions (RHIC) 07J and in relics of 
the early universe and in compact stars |51 Bl 1481 149| . 

Because many theoretical calculations have suggested that the deconfinc- 
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ment transition should be of first order in the low temperature, high density 
regime [SHI OH], we assume it to be a first-order phase transition. If the decon- 
finement transition is of first order, it is natural to believe that in the phase 
transition the system is in a mixed phase. Actually, the hadron-quark mixed 
phase has been considered during hadronization in RHIC (521 031 021 and in 
the transition region from hadron matter to quark matter in neutron stars 

123 123 E3 EH]. 

Our purpose here is to investigate low-density nuclear pasta structure, kaonic 
pasta structure, and hadron-quark pasta structure self-consistently within the 
mean-field approximation. In particular, we figure out how the charge screening 
effects modify the results obtained disregarding these effects in a self-consistent 
manner. 

This review article is organized as follows. In Sec. 2, investigations on the 
mixed phase is reviewed. In Sec. 3 and 4, the nuclear pasta structure in low- 
density nucleon matter and the kaonic pasta structure at higher densities are 
studied. Then the finite-size effects on the pasta structures are discussed in 
Sec. 5. In Sec. 6 and 7, pasta structures in hadron-quark phase transition and 
the charge-screening effects are discussed. Finally, summary and concluding- 
remarks are given in Sec. 8. 

2 General remarks about the treatment of the 
mixed phases 

2.1 Bulk calculation and the finite-size effects 

Consider the mixed phase composed of two different phases denoted by I and 
II. Then the Gibbs conditions require the pressure balance and the equality of 
the chemical potentials between two phases for phase equilibrium |59|. 2 For a 
multi-component system with more than one chemical potential, as is common 
in neutron-star matter, we must impose the equality condition for each chemical 
potential in order to fulfill the condition of the physico-chemical equilibrium. 
More definitely, we hereafter consider the charge chemical potential (/j.q) and 
the baryon- number chemical potential (jj,b) respecting two conservation laws 
in neutron-star matter: /ig = /1q and ^ B — /i^. On the other hand, the first 
condition is not fulfilled in the Maxwell construction, since the local charge 
neutrality is implicitly imposed, while only the global charge neutrality must be 
satisfied. 

A naive application of the Gibbs conditions to the mixed phase composed of 
two semi-infinite matter, when one ignores the surface and Coulomb interaction, 
demonstrates a broad region of the structured mixed phase (SMP), cf. [2"§1 150] . 
However the charge screening effect (caused by the non-uniform charged particle 
distributions) should be very important when the typical structure size is of the 
order of the minimal Debye screening length in the problem. It may largely 

2 We consider here matter at zero temperature 
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affect the stability condition of the geometrical structures in the mixed phases. 

We shall see that the Debye screening effects greatly modify the mechanical 
stability of SMP, and consequently largely limit the density region of the mixed 
phase. In the absence of SMP we effectively recover the picture of phase equi- 
librium given by the Maxwell construction where two bulk phases are separated 
without spoiling the Gibbs conditions. 

Consider SMP consisting of two phases I and II, where we assume spherical 
droplets of phase I with a radius R to be embedded in the matter of phase II 
and two phases are clearly separated by sharp boundaries. We divide the whole 
space into the equivalent Wigner-Seitz cells with a radius i?w (see Fig.|2J). The 
volume of the cell is Vw = iirR^/S and that of the droplet is V = 4ttR 3 /3. 




Vw 



Figure 2: Equal droplets made of the phase I embedded in the phase II (right 
panel), and the geometrical structure of the Wigner-Seitz cell (left panel) 

A bulk calculation proceeds as follows [57]. First, consider two semi-infinite 
matter with a volume fraction fy, separated by a sharp boundary. By applying 
the global charge-neutrality condition and the pressure balance condition under 
the chemical equilibrium, we can get fy for each density. Then we use the 
volume fraction thus determined to describe the geometrical structure in each 
Wigner-Seitz cell. For a given volume fraction fy — (R/Ryj) 3 , the total energy 
E may be written as the sum of the volume energy Ey, the Coulomb energy 
Ec and the surface energy Es, 

E = E V + E C + E S . (1) 

We further assume, for simplicity, that baryon number (pg) and charge (p" h ) 
densities are uniform in each phase a, a = I, II as in semi-infinite matter. 
Then, Ey can be written as Ey /Vw = fv^{p l s) + (1 — /V) eII (Ps) in terms of 
the energy densities e a , a = I, II. The surface energy Es may be represented 
as Es/Vw — fv x 3r/iZ in terms of the surface tension r. The Coulomb energy 
Ec is given by 

W^/.x^^-p^) 2 ^. (2) 

The optimal value of R, which we call i?D> is determined by the minimum 
condition, 

d{E/Vw) 



dR 



= 0, (3) 

R=Ro 
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for a given fy (see Fig. |3J). Since Ey does not depend on R, we can always 
find a minimum as a result of the competition between the Coulomb and the 
surface energies, satisfying the well-known relation, Eg — 2Ec- However, such 




Figure 3: Schematic view of the total energy and each contribution in the bulk 
calculations (solid curves). Screening effect reduces the Coulomb energy, shown 
by the thick arrow. 

bulk calculations have been proved to be too crude for the discussions of SMP. 
Instead, a careful consideration of the interface of two phases is required. As a 
defect of the bulk calculations they ignore the finite-size effects. In particular, 
they have the inconsistent treatment of the Coulomb potential; they do not use 
the Poisson equation, so that the charge density profiles are assumed ab initio 
to be uniform and the Coulomb potential is assumed to be I/r. If one properly 
solves the Poisson equation, one should have the screening effect as a result 
of the rearrangement of the charge-density distribution. Hence, the radius Rd 
should be not too large compared with the Debye screening length Ad, 

V* - % (5) 

in order the above treatment to be justified, the suffix i runs over the particle 
species. Otherwise, the Coulomb energy is reduced by the charge screening ef- 
fect, which should lead to a mechanical instability of SMP in some cases (Fig.[3J). 

We have been recently exploring the effect of the charge screening in the 
context of the various structured mixed phases ED EH ESI Ell E3 ESI- I n 
fact, we have examined the mixed phase at the hadron-quark transition, kaon 
condensation and of nuclear pasta, and found that in cases of the hadron-quark 
transition and kaon condensation the mixed phase might be largely limited 
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by the charge screening and surface effects. In the case of the hadron-quark 
deconfinement transition, for example, A D ~ 5 fm and A D ,A D are of the same 
order as A D , for a typical density with \i& — 1 GeV. We shall see in the following 
that i?D is typically of the same order as Ad ~ A D , and the mechanical stability 
of the droplet is much affected by the screening effect. 

Besides the screening effect, note that the surface tension and the pressure 
balance condition are inconsistently incorporated in the bulk calculations; ac- 
tually the pressure balance condition is imposed before introducing the surface 
tension. 



2.2 Gauge invariance and the meaning of chemical poten- 
tials 

When we use the idea of the density functional theory [S3 ISB]> the thermody- 
namic potential is given as a functional of the particle densities Pi(r), 

Q(pi(r)) = J d 3 re kin+stl .( Pi (r)) + E v -Y^tM J d 3 r Pi (r), (6) 

i 

where pi are the chemical potentials and ekin+str stands for the contributions of 
the kinetic energy and the strong interaction energy. The Coulomb interaction 
energy Ey is also expressed in terms of particle densities, 

e v = -Y ( d 3 rd 3 r '9i£p9i^2 ; (7) 

r - r' w 

with Qi being the particle charge (Q = — e < for the electron). Then the 
chemical potentials are given as 

f M = ^f^-N^V(r), N^ = Q i /e, (8) 

with the electric potential V(r): 

generated by the particle distributions. 

We must keep the gauge invariance through the calculation: V can be arbi- 
trarily shifted by a constant Vq, V(r) — » V(r) — Vq. Formally varying Eq. JHJ) 
with respect to V(r) or pi we have the matrix form relation, 

Ai imr =N ° h > A ^B jk = S ik , (10) 

where matrices A and B are defined as 

B^^i. (ii) 
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Eqs. I|1(J|) and Ijllfl reproduce the gauge-invariance relation, 



dpi 

dv 



N ch9pj 



(12) 



clearly showing that constant-shift of the chemical potential is compensated by 
the gauge transformation of V(r): fa — ► fa+Nf h Vo, as V(r) — > V(r) — Vq. Hence 
chemical potential fa acquires physical meaning only after fixing the gauge of 
V(r). 

We reconsider the relation between the Gibbs conditions and the Maxwell 
construction from this view point. Consider a schematic situation for simplicity 
where two semi-infinite matter denoted by I and II are separated by a sharp 
boundary. As has been mentioned, at first glance the Maxwell construction 
apparently violates the Gibbs conditions, especially the equilibrium condition 
for the charge chemical potential (J.q(= fa) m our context. However, correctly 
speaking, when we say fa\ ^ p]} within the Maxwell construction, it means 
nothing but the difference in the electron number density p e in two phases, 
Pe Pe"> this is simply because p e — pl/(3ir 2 ), if the Coulomb potential is 
absent. Once the Coulomb potential is taken into account, using Eq. (JHJ, p e can 
be written as 

_ (He - Vf 



3tt 2 



(13) 



Thus we may have p\ — pf} and p\ ^ p]} simultaneously, with the different 
values ofV, V 1 ^ V 11 (see Fig.©. We shall see that if the Coulomb interaction 
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Figure 4: Relation between the charge chemical potential P-q{— p- e ) and the 
electron number density p e in the presence of the Coulomb potential V . Fulfilling 
the Gibbs conditions, p) B — /i^, p l e — /ij, 1 , we can change p e in two phases as in 
the Maxwell construction, if V suitably changes from one phase to another. 

is properly taken into account, the resultant EOS looks similar to that given by 
the Maxwell construction. 
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3 Low-density nuclear pasta structures 



3.1 Relativistic mean-field (RMF) treatment of nucleon 
matter 

Here we explain how to investigate the property of nuclear matter at low density. 
Exploiting the idea of the density functional theory within the RMF model, we 
can formulate equations of motion to study non-uniform nuclear matter numeri- 
cally, cf. |^| |S5]. The RMF model with fields of mesons and baryons introduced 
in a Lorentz-invariant way is not only relatively simple for numerical calcula- 
tions, but also sufficiently realistic to reproduce bulk properties of finite nuclei 
as well as the saturation properties of nuclear matter |69) . In our framework, 
the Coulomb interaction is properly included in the equations of motion for 
nucleons and electrons and for meson mean fields, and we solve the Poisson 
equation for the Coulomb potential V self-consistently with those equations. 
Thus the baryon and electron density profiles, as well as the meson mean fields, 
are determined in a fully consistent way with the Coulomb interaction. 



3.1.1 Thermodynamic potential and the equations of motion 



We start with the thermodynamic potential for the system of neutrons, protons, 
electrons and mesons including kaons 



n = n 



N 



M 



The first term 



a— p,n 



(fir 



fcF - <fik_ 

47T 3 



— r-i/TO* +k 2 -p a V a 



(14) 



(15) 



is the contribution of nucleons with the local Fermi momenta fcp,a(r); a = n,p, 
m* N (r) = rnjv — g a Na{v) is the effective nucleon mass and mjy is the nucleoli 
mass in the vacuum. Nucleons couple with a, to and p mesons and thereby, 



v n {r) = fi n - g^AfWo(r) + 5 P jvi?o(r), 

u p (r) = (i p + V(r) - # W ivM)(r) - g pN R (r), 



(16) 



where p n and fi p are neutron and proton chemical potentials and g a N , g^N and 
g P N are coupling constants between corresponding fields. 

The second term in 1)140 incorporates the scalar (a) and vector (u>o, Rq) mean 
fields, 



n 



M 



d 6 r 



(Vcr) 2 + mla 2 



(Vcjq) 



Uia) 



2 < mlujl (Vi? ) 2 + mlRl 



(17) 



10 



Table 1: Parameter set used in RMF in our calculation. 



9crN 


9ujN 


9pN 


b c m a m u m p 


6.3935 


8.7207 


4.2696 


0.008659 0.002421 400 MeV 783 MeV 769 MeV 



where m CT , m w and m p are the field masses, and U(a) = ^bm^igaN^) 3 + 
\c(g a N<j) A is the nonlinear potential for the scalar field. 

The third term in (|14|> contains the contribution of the Coulomb field (de- 
scribed by the potential V(r)) and the contribution of relativistic electrons, 



-;3, 



1 -_ ma fa-V) 



;(W) 



8ire 2 x ' 12tt 2 



(18) 



where u, e is the electron chemical potential. 

Temperature T is kept to be zero in the present study. 

For nucleons and electrons we used the local-density approximation, i.e., nu- 
cleons and electrons are described by their local densities. This approximation 
has its sense only if the typical length of the change of the nucleon density is 
larger than the inter-nucleon distance. Derivative terms of the particle densities 
can be incorporated in the quasi-classical manner by the derivative expansion 
within the density functional theory jB3|HS]. Their contribution to the energy 
can be reduced to a surface tension term. Here we simply discard those deriva- 
tive terms, as a first-step approximation. Thus we discard the contribution of 
the nucleon fields to the surface tension assuming that it is smaller than the 
corresponding contribution of the meson fields that we retain. In the case when 
we suppress derivative terms of the nucleon densities they follow changes of the 
meson a, lu, p mean fields and the Coulomb fields that have derivative terms. 
Note that we have fitted our model to properly describe finite nuclei (see below) 
without including nucleon derivative terms. If we introduced them it would 
need to re-adjust the model parameters. We should bear in mind that for small 
structure sizes quantum effects become prominent. For simplicity we disregard 
these effects. Thus we may properly describe only rather large-size structures 
within this scheme. 



3.1.2 Equations of motion and numerical procedure 

Equations of motion for the mean fields and for the Coulomb potential are 

so. 



obtained from the variational principle: , f 9 \ = (<j>i = cr,ujo, Ro, V), 



V 2 a(r) = mJaW + l-j^rJ+^r)), (19) 
V 2 w (r) = r7i 2 wo(r) -g U N(p P (r) +P„(r)), (20) 
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V 2 i? (r) - m 2 p R (r)~g pN (p p (r)~ p n (v)), (21) 
V 2 V{v) = 4^e 2 P ch(r), (22) 

with the proton and neutron scalar densities 

f kF ' a(r) d3k m * N (r) 
Jo 47r V m Jv( r ) + k 



and the charge density 



Pch(r) = ft,(r)-p e (r). (24) 



For nucleons and electrons, Sp ^ — (a = n,p,e) give the expressions of the 
chemical potentials; 

Me = (37r 2 p e (r)) 1/3 + y(r), (25) 

Mn = MB = \/^F,n(r) 2 +TO^(r) 2 

+ 5wArWo(r) - g pN R (r), (26) 

M P = Mb - Me = \/ fcF,p(r) 2 + m^(r) 2 

+ 3a ) JV^'o(r) + 9 P NRo(r) ~ V(r). (27) 

The last two equations are the standard relations between the local nucleon 
densities and chemical potentials. We have assumed that the system is in the 
chemical equilibrium with respect to the weak, electromagnetic and strong in- 
teractions and we introduced the baryon chemical potential Ms = M« an d the 
charge chemical potential, i.e. the electron chemical potential, /i e , according to 
the corresponding conserved charges. 

To solve the above coupled equations numerically, the whole space is divided 
into equivalent Wigner-Seitz cells with a radius The geometrical shape of 
the cell changes as follows: sphere in three-dimensional (3D) calculation, cylin- 
der in 2D and slab in ID, respectively. Each cell is globally charge-neutral and 
all physical quantities in the cell are smoothly connected to those of the next cell 
with zero gradients at the boundary. Every point inside the cell is represented 
by the grid points (number of grid points -/V gr id w 100) and differential equations 
for fields are solved by the relaxation method for a given baryon-number den- 
sity under constraints of the global charge neutrality. Details of the numerical 
procedure are explained in Ref. (HHI- To fix the gauge of V^r), we choose the 
following condition as 

V(R W ) = 0. (28) 



3.1.3 Parameter set and finite nuclei 

Parameters of the RMF model are chosen to reproduce saturation properties 
of symmetric nuclear matter: the minimum energy per nucleon —16.3 MeV at 
p = po = 0.153 fm~ 3 , the incompressibility K(po) = 240 MeV, the effective 



12 




Figure 5: EOS of uniform nucleon matter with different proton mixing ratios. 
Red line indicates symmetric matter (proton mixing ratio Y p — 0.5) and blue 
line indicates neutron matter (Y p = 0). Lines with intermediate colors show 
the case of Y p — 0.4, 0.3, 0.2 and 0.1. Crosses on the lines show the minimum 
points. Upper panel is the case without electron and the Coulomb energy, while 
the lower panel is the case of charge-neutral matter which includes energy of 
electrons. 

nucleon mass m* N {po) — 0.78toat; tun — 938 MeV, and the isospin-asymmetry 
coefficient a sym = 32.5 MeV. Coupling constants and meson masses used in 
our calculation are listed in Table 1. Figure shows the binding energy per 
nucleon of uniform nucleon matter. The proton mixing ratio Y p is the number 
ratio of protons (Z) to total baryons (A), Y p = Z/A. The saturation property 
mentioned above is seen for symmetric nuclear matter (Y p — 0.5) in the upper 
panel of Fig. Note that if one imposes the charge neutrality by inclusion of 
electrons, the saturation property cannot be seen (lower panel of Fig.[5J). 

For the study of non-uniform nuclear matter, the ability to reproduce the 
bulk properties of finite nuclei should be essential. We check how it works to 
describe finite nuclei. In this calculation, for simplicity, we assume the spherical 
shape of nuclei. The electron density is set to be zero. Therefore, neither the 
global charge neutrality condition nor the local charge-neutrality condition is 
imposed. 

In Fig. [5] (left panel) we show the density profiles of some typical nuclei. One 
can see how well our framework may reproduce the density profiles. To get a 
still better fit, especially around the surface region, we might need to include 
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Figure 6: Properties of finite nuclei obtained with the present RMF model. Left: 
the density profiles of typical nuclei. The proton number densities (solid curves) 
are compared with the experimental data. Right: binding energy per nucleon 
— (E/A — tun) and the proton mixing ratio Y p of finite nuclei. 

the derivative terms of the nucleon densities, as we have already remarked. Fine 
structures seen in the empirical density profiles, which may come from the shell 
effects (see, e.g., a proton density dip at the center of a light 16 nucleus), cannot 
be reproduced by the mean-field theory. The effect of the rearrangement of the 
proton density distribution is seen in heavy nuclei. Protons repel each other, 
which enhances their density near the surface of heavy nuclei. This effect is 
analogous to the charge screening effect for the Coulomb potential in a sense that 
the proton distribution is now changed not on the scale of the nuclear radius, but 
on another length scale, that we will call the proton Debye screening length (see 
Eq. 130|) below). It gives rise to important consequences for the pasta structures 
since the proton Debye screening length is typically less than the lump size. The 
optimal value of the proton mixing ratio Y p is obtained by imposing the beta- 
equilibrium condition for a given baryon number. Figure (right panel) shows 
the baryon-number dependence of the binding energy per baryon -Ebind and the 
proton mixing ratio Y p . We can see that the bulk properties of finite nuclei 
(density, binding energy, and proton to baryon number ratio) are satisfactorily 
reproduced for our present purpose. 

Note that in our framework we must use a sigma mass m a = 400 MeV [20] > 
a slightly smaller value than that one usually uses, to get an appropriate fit. 
If we used a popular value m a w 500 MeV, finite nuclei would be over-bound 
by about 3 MeV/ A. The actual value of the sigma mass (as well as the omega 
and rho masses) has little relevance for the case of infinite nucleon matter, since 
it enters the thermodynamic potential only in the combination C a = g a N /m a . 
However meson masses are important characteristics of finite nuclei and of other 
non-uniform nucleon systems, like those in pasta. The effective meson mass 
characterizes the typical scale for the spatial change of the meson field and 
consequently it affects the value of the effective surface tension. 
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Figure 7: Examples of the density profiles in the cell for symmetric nuclear 
matter with 1^=0.5 (left panel) and for asymmetric matter with Y p = 0.3 (center 
panel) and 0.1 (right panel). 

3.2 Nucleon matter at fixed proton mixing ratios 

First, we concentrate on the discussion of the behavior of nucleon matter at 
fixed values of the proton mixing ratio Y p . Particularly, we explore the proton 
mixing ratios Y p — 0.1, 0.3, and 0.5. The cases Y p = 0.3 - 0.5 should be 
relevant for supernova matter and for newly born neutron stars. Figure[3shows 
some typical density profiles inside the Wigner-Seitz cells. The geometrical 
dimension of the cell is denoted as "3D" (three-dimensional sphere), etc. The 
horizontal axis in each panel denotes the radial distance from the center of the 
cell. The cell boundary is indicated by the hatch. From the top to the bottom 
the configuration corresponds to droplet (3D), rod (2D), slab (ID), tube (2D), 
and bubble (3D). The nuclear "pasta" structures are clearly manifested. For 
the lowest Y p case (Y p = 0.1), the neutron density is finite at any point: the 
space is filled by dripped neutrons. The value of Y p above which neutrons drip 
is around 0.26 in our 3D calculation, for example. For a higher Y pi the neutron 
density drops to zero outside the nucleus. The proton number density always 
drops to zero outside the nucleus. We can see that the charge screening effects 
are pronounced. Due to the spatial rearrangement of electrons the electron 
density profile becomes no more uniform. This non-uniformity of the electron 
distribution is more pronounced for a higher Y p and a higher density. Protons 
repel each other. Thereby the proton density profile substantially deviates from 
the step-function. The proton number is enhanced near the surface of the lump. 

The equation of state (EOS) for the sequence of geometric structures is 
shown in Fig. |H1 (top panels) as a function of the averaged baryon-numbcr 
density. Note that the energy E/A — uin also includes the kinetic energy of 
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Figure 8: Binding energy per nucleon (top), pressure (middle), and the cell and 
lump sizes (bottom) for symmetric nuclear matter with Y p =0.5 (left panel), and 
for asymmetric matter with Y p = 0.3 (center panel) and 0.1 (right panel). 



electrons, which makes the total pressure positive (middle panels). The lowest- 
energy configurations are selected among various geometrical structures. The 
most favorable configuration changes from the droplet to rod, slab, tube, bub- 
ble, and to the uniform one (the dotted thin curve) with an increase of density. 
The appearance of non-uniform structures in matter results in a softening of 
EOS: the energy per baryon gets lower up to about 15 MeV/^4 compared to 
uniform matter. 

The bottom panels in Fig. [SI are the cell radii i?w and the lump radii R 
versus averaged baryon number density. The radius R is defined by way of a 
density fluctuation as 

Rw P 2 , i (for droplet, rod, and slab) 

/> = { f p) (p)2 , (29) 
Rw 1 - t^tt ' ( for tube and bubble ) 



(P 2 P ) 

where the bracket "()" indicates the average along the radial (for 3D and 2D 
cases) or perpendicular (ID) direction in the cell. Dashed curves show the Debye 
screening lengths of electrons and protons calculated as 



A ° = ( a =P> e )> ( 3 °) 

where Pp V is the proton number density averaged inside the nucleus (the region 
with finite p p ) and pf is the electron density averaged inside the cell. Actually 



-1/2 



1G 



doing more carefully we should introduce four Debyc screening lengths A^ :< and 
with a separate averaging for the interior and the exterior of the nuclei. 
However we observe that the proton number density is always zero in the exte- 
rior region and =00 thereby. For electrons A^ < and A^ are in general 
different but both being large and of the same order of magnitude in the pasta 
case under consideration. Therefore we actually do not need a more detailed 
analysis of these quantities. Note that these values are obviously gauge invari- 
ant. Numerically the cell radii i?w for droplet, rod, and slab configurations at 
Y p = 0.5 and 0.3 were proven to be close to the electron screening length. For 
the tube, i?w is larger than Ap. For Y p = 0.1, in all cases i?w is substantially 
smaller than Xfy and the electron screening should be much weaker, thereby. In 
all cases, except for bubbles (at Y p = 0.5 and 0.3), the structure radii R are 
smaller than A£> . This means that the Debye screening effect of electrons inside 
these structures should not be pronounced. For bubbles at Y p = 0.5 and 0.3, A^ 
is substantially smaller than the cell size and the electron screening should be 
significant. For Y p — 0.5, 0.3, 0.1 in all cases (with the only exception Y p = 0.1 
for slabs), the value Aj^ is shorter than R. Hence the density rearrangement of 
protons is essential for the pasta structures, as it is indeed seen from the Fig.0 
As we have mentioned one of the characteristics of our calculation is that 
the obtained density profiles of particles in the cell is free from assumption and 
is consistent with the potential. However, the usage of Wigner-Seitz cell is a sig- 
nificant limitation. It is because we have to assume the geometrical symmetry. 
At zero temperature, a specific pasta structure should be the distinct ground 
state for a given baryon density. However, at a certain density near the bound- 
ary of two pasta structures, some intermediate structures could be favorable, 
or at least can be local-minimum state. Furthermore, intermediate structure 
should not be negligible at finite temperatures. Such intermediate states and 
incomplete pasta structures are reported in the studies using molecular dynam- 
ics ^3 1201 Ell E2 • In fact, molecular dynamics simulation is completely free 
from assumption of nucleoli distribution. However, uniform electron background 
is assumed and beta-equilibrium condition is not imposed in the molecular dy- 
namics yet. In this sense, the present mean-field calculation and the molecular 
dynamics for nuclear matter structure are complementary models to each other. 

3.3 Nucleon matter in beta equilibrium 

Next, we consider neutron star matter at zero temperature, and explore the 
non- uniform structures for nucleon matter in beta equilibrium. Figure El shows 
the density profiles for different baryon number densities. The droplet structure 
itself is quite similar to the case of the fixed proton mixing ratio Y p = 0.1 
considered above. The apparently different feature in this case is that only the 
droplet configuration appears as a non-uniform structure. It should be noticed, 
however, that the presence or absence of the concrete pasta structure sensitively 
depends on the choice of the effective interaction. 

In Fig. we plot the energy per baryon (top), the cell and lump sizes 
(middle), and the proton number ratio (bottom). The effect of the non-uniform 
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Figure 9: Density profiles in the cell 
for nuclear matter in beta equilibrium 
with baryon-number densities, 0.01, 
0.03 and 0.05 fm -3 from the top to 
the bottom. 




Figure 10: Binding energy (top), cell 
and lump sizes (middle), and proton 
mixing ratio (bottom) in the cell for 
nuclear matter in beta equilibrium. 



structure on EOS (the difference between the energy of uniform matter and that 
of non-uniform one) is small. However, the proton mixing ratio is significantly 
affected by the presence of the pasta at lower densities. In the zero-density 
limit, the proton mixing ratio should converge to that of the normal nuclei. 
The droplet radius and the cell radius in the middle panel of Fig. E3 are always 
smaller than the electron Debye screening length Ap . Thereby the effect of the 
electron charge screening is small. Unlike the fixed Y p cases, the droplet radius 
is comparable to the proton Debye screening length, which means that the effect 
of the proton rearrangement is not pronounced in this case. In fact, there is 
no enhancement of the proton number density near the surface in Fig. [5J in 
contrast to Fig. 



4 Kaonic pasta structures at high densities 

Next let us explore high-density nuclear matter in beta equilibrium, which is 
expected in the inner core of neutron stars. 

Kaons are Nambu-Goldstone bosons accompanying the spontaneous break- 
ing of chiral SU(3) x St/ (3) symmetry and the lightest mesons with strangeness. 
Their effective energy is much reduced by the kaon-nucleon interaction in nu- 
clear medium, which is dictated by chiral symmetry. For low-energy kaons the 
s-wave interaction is dominant and attractive in the I — 1 channel, so that 
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Figure 11: Onset mechanism of kaon 
condensation within the chiral model. 



Figure 12: Original EOS and the 
Maxwell construction given by the chi- 
ral model. "C" indicates the criti- 
cal point between nucleon and kaonic 
phases in the case of uniform matter 
(without mixed phase). The mixed 
phase is given by the Maxwell con- 
struction between "N" and "K" . 



negatively charged kaons appear in the neutron-rich matter once the process 
n — > p + K~ becomes energetically allowed. Since kaons are bosons, they cause 
the Bose-Einstein condensation at zero momentum j21j . The single-particle 
energy of kaons is given in a model- independent way (see Appendix) : 

e±(p) = V / |p| 2 + mf + ((p„ + 2 Pp )/4/2)2± ( Pn + 2p p )/lf, (31) 

where m* K is the effective mass of kaons, m*^ = m 2 K — ^kn(p u + p s p )/ f 2 , with 
the KN sigma term, E^tv, and the meson decay constant, / = fx ~ f n . The 
threshold condition then reads [73] 

p K = e_(p = 0) = p n - p p = p e , (32) 

which means the kaon distribution function diverges at p = fFig. Hip . 

If kaon condensation occurs in nuclear matter, it has many implications on 
compact stars; one of the interesting possibilities may be the delayed collapse of 
a protoneutron star to a low-mass black hole due to the large softening of EOS 
|32l I35| . and another one is a fast cooling mechanism of young neutron stars 
due to the nucleon Urea process under background kaons |361 137| . 

Since many studies have shown that kaon condensation is of first order, we 
must carefully treat the phase change fFig. 1121) . In the following we discuss the 
mixed phase in the phase transition in a similar way to nuclear pasta. 
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4.1 RMF treatment of nuclear matter with kaon conden- 
sation 



We explore high-density nuclear matter with kaon condensation by means of 
RMF model as in low-density matter. Using the same model we can discuss the 
non-uniform structure of nuclear matter both at low- and high-density regime 
in an unified way. To incorporate kaons into our RMF calculation, the thermo- 
dynamic potential of Eq. 114fl is modified as 



n 



n 



K 



where 



N 



d 3 r 



M 

f 2 K o 2 r 



K ! 



' K 



(33) 



+ (fi K - V + g U KUo + g P KRoY 



nm 2 



(34) 



'K 



rriK - gaKcr, hk 



e , and the kaon field K = JkO/V^ (Jk is 
the kaon decay constant). We, hereafter, neglect a rather unimportant term 
oc a 2 9 2 . 3 The equations of motion are similar to Eqs. I|19|) - l|27|> given for the 
low-density case (Sec. 3). We list here the modified and added ones as 



V 2 a(r) 

V 2 w (r) 
V 2 i? (r) 



mla{r) + ^- - g<r N (p s n (r) + P 6 P ( r )) 

+ 2g aK m K f K e 2 (r), 
mlu) (r) ~ g wN {p P {*) + Pn(r)) 



(35) 



V 2 6»(r) 



- f K 9ujKd 2 (r)[fi K - V(r) + g uK Uo(r) + g pK Ro(r)}, (36) 
m 2 p R (r) - g pN (p p {v) - p n (r)) 

- f 2 K9 P K0 2 {Y)[pK - V(r) + g uK uo(r) + g pK R {r)}, (37) 
m* K (r) 2 - {\i K - V(r) 

+ 9u K u (r) + g pK R Q {r)) 2 } 9(r). (38) 
The charge density now reads 

Pch(r) = Pp(r) ~ Pe(r) - Asr(r), (39) 
with the kaon contribution, 



Pk = f K 8 [fJ-K -V + g uK u (r)+g pK R (r)} 



(40) 



Additional parameters concerning kaons are presented in Table 2. By a lattice 
QCD calculation, the kaon-nucleon sigma term T,kn is estimated as 290 - 450 
MeV [7*!]. Using the above values and the relation (see Appendix) 



g<rKg<rN 



2m K fK 



2 ' 



3 We use the meson exchange model (MEM) for kaon Lagrangian to compare our results 
with the earlier ones 28 30 31 , which is not chiral-symmetric (see Appendix). Note however 
that both can well describe the low-energy kaon-nucleon dynamics. 
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Table 2: Additional parameters used in our RMF model with kaon terms. The 
kaon optical potential Uk is defined by Uk = goKO + JuKWo- 

U) [MeV] m K [MeV] ~g~K Efc(po) [MeV] 



93 494 ff^jy/3 g pN -120 - -130 



9<tK can be 0.849 - 1.318. On the other hand, the parameter g a K enters the 
value of the K~ optical potential Uk defined by Uk = daK® + g^K^o- There 
have been many works trying to extract Uk at the saturation density from the 
data on the kaonic atoms [75] [7Sj and from calculations %77\ 1781 1791 1801 18 1| . 
but there is still a controversy about its depth. We take here a somewhat 
deep potential, as shown in Table 2, to compare our results with the earlier 
ones HD ED ED EI]- The corresponding g aK can be 2.209 - 2.519, which may 
be rather large compared to the lattice QCD estimation. To understand a 
dependence of the results on the value of Uk we further allow for its variation. 

4.2 Properties of kaonic pasta structures 

The kaon condensation is considered to be of first order. Therefore it may give 
rise to non-uniform structures with different geometries, the structured mixed 
phase. In fact, the system exhibits a series of structure change similar to the nu- 
clear "pasta": the kaonic droplet, rod, slab, tube, bubble (which we call "kaonic 
pasta" structures). Figure ITTA displays typical density profiles and the Coulomb 
potential. The neutron distribution proves to be rather flat. The proton dis- 
tribution on the other hand is strongly correlated with the kaon distribution, 
which means that the Coulomb interaction is crucial. 

In the upper panel of Fig. 1141 we depict the energy per nucleon as a function 
of baryon number density. The dotted line indicates the case of single phase 
(if one assumes the absence of the mixed phase). In this case uniform matter 
consists of normal nuclear matter below the critical density and kaonic matter 
above the critical density. The cross on the dotted line (ps — 0.46 fm -3 ) 
shows the critical density, i.e. the point where kaons begin to condensate in 
the case of single phase. Pieces of solid curves, on the other hand, indicate 
energetically favored structures. Droplets begin to appear for ps > 0.41 fm -3 
and smoothly decrease the energy of the system. The mixed phase disappears 
for p B > 0.74 fm~ 3 . 

The lower panel of Fig. 1141 shows the sizes of the lump R and the cell i?w- 
We find that at the onset density the size of the cell is infinitely large in case 
of the full calculation. The corresponding steep increase of i?w with decreasing 
density is clearly seen in the figure. The dashed lines and the dotted line in the 
lower panel of Fig.^Jshow partial contributions to the Debye screening lengths 
of the electron, proton and kaon, Aq, Aq, and Ap , respectively. We see that in 



21 



l [3D, pg =0.45fm" 3 
0.5 



Ug=-130 MeV 




Figure 13: Density profiles of kaonic 
pasta structures. Here the density 
does not mean charge-density but 
number-density of particles. Uk, the 
kaon optical potential at the nuclear 
saturation density, is set to be —130 
MeV. Displayed using the right axis is 
the Coulomb potential V. 



Figure 14: Top: binding energy per 
nucleon of the nuclear matter in beta 
equilibrium. The dotted line below the 
cross shows uniform normal nuclear 
matter and above the cross, uniform 
kaonic matter. Bottom: lump size R 
(thick curves below) and cell size i?w 
(thick curves above). Compared are 
the Debye screening lengths of elec- 
tron, proton and kaon. The Debye 
screening lengths are calculated using 
the explicit dependence of fi a on p a 
(see Eq. (GP). 
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Figure 15: Top: binding energy per nucleon of nuclear matter. Bottom: lump 
size R and cell size i?w- The proton mixing ratio is Y p =0.5 for all cases. Different 
calculations are compared. From the left: "full", "no Coulomb" and "weak 
surface" calculations. In "no Coulomb" calculation, the electric potential is 
discarded when determining the density profile and then added to evaluate the 
energy. 

most cases is less than the cell size i?w but is larger than the lump size R. 
The proton Debye length and the kaon Debye length A^ , on the other hand, 
are always shorter than Rw and R. 
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Figure 16: The same as Fig. 1151 for high-density nuclear matter in beta equi- 
librium. 
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5 Charge screening and surface effects on the 
pasta structures 

5.1 Charge screening effect 

To demonstrate the charge screening effects we compare results of the full calcu- 
lation with those given by a perturbative treatment of the Coulomb interaction 
often used in the literature, which we refer to as "no Coulomb" calculation. In 
this calculation the electric potential is discarded in equations of motion (|25|) , 
JSJ, (jlUjl which determine the density profiles. The Coulomb en- 

ergy is then added to the total energy by using the charge density profile thus 
determined to find the optimal value with respect to the cell size i?w 4 

First, we discuss the case of low-density symmetric nuclear matter. In the 
left and the central panels of Fig. compared are different treatments of 
the Coulomb interaction. The EOS (upper panels) as a whole shows almost no 
dependence on the treatments of the Coulomb interaction. However, sizes of the 
cell and the lump (lower panels), especially for tube and bubbles, are different. 
In the cases of the "full calculation", the cell radii of "tube" and "bubble" 
structures and that of "slab" structure get larger with increase of density, while 
they are monotonically decreasing in the case of "no Coulomb" calculation. The 
other effect is a difference in the density range for each pasta structure. The 
"full" treatment of the Coulomb interaction slightly increases the density region 
of the nuclear pasta. 

We show the same comparison for the kaonic pasta structures in Fig. 1161 We 
see again that the density range of the mixed phase is narrower in the case of 
the "no Coulomb" calculation than in the full calculation, while the energy gain 
is almost the same. A remarkable difference is seen in the cell size, especially 
near the onset density of kaonic pasta, for pb < 0.5 fm -3 . The cell size given 
by the full calculation is always larger than that given by the "no Coulomb" 
calculation. 

To elucidate the screening effect, we depict the Rw dependence of the energy 
per nuclcon in Fig. El ln a general case of 3D droplet the Coulomb energy 
per particle depends on the radius by its square, while the surface energy per 
particle by its inverse. Therefore the sum of the Coulomb and surface energy 
has a U-shape (cf. "no Coulomb") and has a minimum at a certain radius. If the 
Coulomb interaction is screened, the Coulomb contribution will be suppressed 
(cf. the full calculation) and the minimum point gets larger (see Fig. |3J). Since 
the cell radius is approximately proportional to the droplet radius for a given 
baryon density, the above argument applies also to the cell size. 

For a long time there existed a naive view that not all the Gibbs conditions 
can be satisfied in a description by the Maxwell construction if there are two or 
more independent chemical components [SHUSH EDI) because the local charge 

4 This treatment may look similar to the bulk calculation. However, the pressure balance 
condition is correctly imposed in this calculation, while it is inconsistently used in the bulk 
calculation. 
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Figure 17: The cell size Rw dependence of the energy per nucleon. Crosses on 
the curves indicate the minimum points. 

neutrality is implicitly assumed in it. As the result of this argument, it was sug- 
gested that a broad region of the structured mixed phase may appear in neutron 
stars. However, in recent papers j^E21ESlESlEl we have demonstrated that if 
one properly includes the Coulomb interaction, the Maxwell construction prac- 
tically satisfies the Gibbs conditions and the range of the mixed phase will be 
limited. 

We present in Fig. I18l the pressure as a function of the baryon-number den- 
sity. We also depict the pressure when the Gibbs conditions are applied for two 
semi-infinite matters discarding the Coulomb interaction (indicated by "Gibbs" ) 
and that given by the Maxwell construction (indicated by "Maxwell" ) . We see 
that the pieces of solid curves lie between "Gibbs" and "Maxwell" . The full cal- 
culation case is more similar to the one given by the Maxwell construction. In 
Fig. El compared are the density profiles obtained by the full and "no Coulomb" 
calculations. In case of the full calculation the difference between the negative 
charge density (of kaons and electrons) and the positive charge density of protons 
is smaller, indicating that the system tends to have a local charge neutrality. 
These results suggest that the Maxwell construction is effectively justified in the 
full calculation owing to the charge screening effects. 

5.2 Surface effects 

If we artificially multiply meson masses m a , m u and m p by a factor cm, e.g. 
cm = 1 (realistic case), 2.5 and 5.0, the surface tension changes. Figure EH 
demonstrates the binding energies of finite nuclei calculated with different meson 
masses. By the use of heavy meson masses, the binding energy of finite nuclei 
(for finite A) approaches to that of nuclear matter indicated by a thick gray 
line. This shows that the surface tension is reduced with increase of the meson 
masses, cf. jH2]- Notice that this statement is correct only if we fix the ratio 



Using the above modified meson masses, we explore the effects of surface 
tension in the following. Right panels of Figs . ITSI and ITol are EOS and the sizes 
of the cell and the lump, but now for the case of an artificially suppressed surface 
tension (cm = 5.0). Comparing them with the left panels, we see that there 
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Figure 20: Mass number A dependence of the binding energy of finite nuclei 
-E-bind = —(E/ A — ttln) without the Coulomb interaction. Thick gray line shows 
the case of nuclear matter. 
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is almost no difference in the EOS. However, there are two differences in the 
case of low-density nuclear matter. First, the density range of pasta structure 
is slightly broader for weaker surface tension. Secondly, the cell size with a 
normal surface tension is larger than the case of weaker one. It means that 
weaker surface tension and stronger Coulomb repulsion cause the similar effects 
on the cell size since the pasta structure is realized by the balance of the both. 

In the case of kaonic pasta, the meson masses have very small effects. The 
a, u> and p mesons have less contribution to the surface tension of kaonic pasta 
but K-N interaction is dominant. 

6 Quark pasta structures in the hadron-quark 
mixed phase 

Hereafter, we consider the phase transition from nuclear matter to three-flavor 
quark matter in beta equilibrium at zero temperature. Accordingly, we must in- 
troduce many chemical potentials for the particle species, but the independent 
ones in this case are reduced to two, the baryon-number chemical potential 
/xb and the charge chemical potential due to the beta equilibrium and the 
condition of total charge neutrality. We can easily see that they are equal to 
the neutron and electron chemical potentials, [i n and /z e , respectively. In the 
mixed phase, these chemical potentials should be spatially constant to satisfy 
the Gibbs condition (GC), while the density profiles arc non-uniform. When 
we naively apply the Maxwell construction (MC) to get the EOS in thermo- 
dynamic equilibrium, we immediately notice that is constant in the mixed 
phase, while \i e is different in each phase because of the difference in the elec- 
tron number of these phases. This is because the MC uses the EOS of bulk 
matter in each phase, which is locally charge-neutral, uniform matter; many 
electrons are needed in hadron matter to compensate for the positive charge of 
the protons, while in quark matter, total charge neutrality is almost satisfied 
without electrons. We shall see how this situation is changed in the presence 
of the Coulomb interaction. Actually we shall see that MC can be effectively 
used, while it apparently violates the Gibbs conditions. 

The treatment of the mixed phase following the dcconfincmcnt transition 
is somewhat different from that in previous sections due to the property of 
the vacuum. When we consider the mixed phase, hadrons reside in the non- 
perturbative vacuum, while quarks in the perturbative one. So the change of 
the vacuum should be also taken into account in this case, which should greatly 
affect the properties of the interface between the hadron and quark phases. 
Unfortunately, since dynamics of the deconfincmcnt transition has not been 
well understood yet, especially in the finite-baryon number density case, we 
must have recourse to some models. We use the MIT bag-model picture in 
our study, where the non-perturbative and perturbative vacua are separated by 
a sharp boundary and the surface is endowed with the surface tension. Note 
that surface tension is not needed in a consistent calculation, because the non- 
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uniform density profiles naturally give its effect. 
6.1 Thermodynamic potential 

We consider the structured mixed phase (SMP) in which one phase is embedded 
in the other phase with a certain geometrical form. As in the previous studies of 
low-density nucleon matter and high-density kaonic matter, we use the Wigner- 
Seitz cell with a size i?w- A cell includes an embedded quark phase with a size 
R. 

The quark phase consists of u, d, and s quarks and electrons. The hadron 
phase consists of protons, neutrons and electrons. We incorporate the MIT bag 
model and assume a sharp boundary at the hadron-quark interface. We use 
the idea of density functional theory (DFT) and incorporate the local density 
approximation (LDA) |67l lf>8| . 

We consider the total thermodynamic potential: 

^total = ^H + ^Q + ^r, (41) 

where J7h and J7q represent the contributions of the hadron and quark phases, 
respectively. We here introduce the surface contribution fl T in order to prop- 
erly treat the geometrical structure. It may be connected with the confining 
mechanism, but unfortunately we have no definite idea to incorporate it. There 
have been many estimations about the surface tension in lattice QCD [831 184) . 
in shell-model calculations 1861 157) and in model calculations based on the 
dual-Ginzburg Landau theory |88) . Here, we parametrize the surface tension 
by r. Then we can write J7 T = tS, with S being the area of the interface. 
Moreover, there may also be a contribution from the curvature term |86| . and 
we regard its effect to be renormalized into r for simplicity. Actually, many 
authors have treated the strength of r as a free parameter and determined how 
the results depend on its value [SHI E3 EFJ • Here we also employ this approach 
by introducing the surface tension parameter r to simulate the surface effect. 

We now present the expressions of thermodynamic potentials. First, the 
Coulomb interaction energy is expressed in terms of particle densities, as 

Ev = \Y.I dW QiP f )Q ^ (r '\ (42) 

where i, j — u,d, s, p,n,e, with Q$ being the particle charge (Q = — e < for the 
electron) and Vvv the volume of a Wigner-Seitz cell. Accordingly, the Coulomb 
potential is defined as 

f(r)=-£/ +V , (43) 

• JVw l r ~ r I 

where Vo is an arbitrary constant representing the gauge degree of freedom. We 
here fix the gauge by stipulating the condition V(Rw) — as before (see Sec. 
2.2). 
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From the above considerations, we obtain the electron contribution and the 
Coulomb interaction energy (in both phases) as 



d 3 r 

Vw=Vq+Vh 



(W(r)) 2 + e e ( Pe (r)) - MeP e (r) + V(r)p e (r) 



= ^ m + nf m , (44) 

4 

(37T 2 p e (r))^ 

where e e (p e (r)) = ^ — ° $ ' is the kinetic energy density of electrons. The 
contributions in the two phases, f2^n anc ^ are obtained by simply dividing 
!] cm by the volumes Vq and Vh, respectively. 

Secondly, in the quark phase, the u and d quarks are treated as massless 
particles, and only the s quark is massive, with m s = 150 MeV. The kinetic 
energy density of the quark of flavor / is simply expressed as jSU] 

3 

e/kin = -^ m j [ x fVf (2z/ + 1) - In (x f + r)f)] , (45) 

where m/ is the quark mass, Xf — pF/(r)/w/, with the Fermi momentum 

p F/ (r) = (7r 2 p/(r))£, and 77/ = ^Jl + xj. 

For the interaction energy, we take into account the leading-order contribu- 
tion coming from the one-gluon exchange interaction: 

e/Fock = ~-^ m j 1 4 - \ [ x fVf ~ ln 0/ + Vf)]j • (46) 

Including this interaction, the quark contribution to the thermodynamic poten- 
tial is expressed as 



n 



Q 



a u + n d + n s + [ d^B + n^, (47) 

JVq 

n f = [ <Pr[e,(pf(r))-n f pf(T)-N,V(T)p f (T)], N f = & (48) 
Jv Q e 

where the energy density e/(p/(r)) represents e/ki n + e/Fock of the / quark, and 
B is the bag constant. The bag constant is taken as 120 MeV/fm 3 , and the 
QCD fine structure constant as a c = 0.4, which were also used by Heiselberg et 
al. [SZj and in Refs. pH RTTl l6l] . 

Thirdly, we consider the hadron contribution. The thermodynamic potential 
for the non-relativistic nucleons becomes 



n H = E N - V fi a f d 3 rp a {v)- f d 3 rV(r) Pp (r) + nf m , (49) 

a=p,n J V" 



where E-^ is the energy of the nucleons, 



£ N = / d 3 r 



£ T^-(3n 2 y pl(r)+e pot (p P (r),p n (r)) 



10m 

a—p,n 



(50) 
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Here we use the effective potentiai e po t(/9 p (r), p„(r)) parametrized by the local 
densities for simplicity, 



tpot 



(r) 



So 
K 



(p n (r) - p p (r)Y 
Po(r) 

(Pn(r) +p P (r)) ( p n (r)+p p 
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(pn(r) + p P (r)) ebind 

2 



1 



PO 



+ C sat (p n (r) + Pp (r)) ( > " (r)+Pp(r) - I 



(51) 



where So, Kq, ebind and C sa t are adjustable parameters which are chosen to 
reproduce the saturation properties of nuclear matter [HOI EH E3- We use 
So = 18 MeV, Ko = 285 MeV, and C sa t = -14.3 MeV at the saturation density 
po- The value ebind — ebind + - 2p " 1 e p (ps = poj^ = 1/2) provides the 
empirical binding energy ebind = —15.6 MeV 

We consider chemical equilibrium at the hadron-quark interface as well as 
in each phase, in which case we have the following: 



Pu + Pe = Pd, Pd = Psi 

Pp + Pe = pn = PB i Pn = pu 

p p = 2p u + pd. 



^Pd, 



(52) 



6.2 Numerical procedure 

We get the equations of motion from the relation l5 ^ tal = (<f>i — p u ( r ), Pd{ Y )i Ps( r ), 
p p (r), p n (r), p e (r), V(r)). The Poisson equation then reads 

V 2 V{y) = 4^e 2 P ch(r) (53) 

where the charge density p^(r) is given by p ch (r) = §p u (r)- gPd(r)- ±p s -p e (r) 
in the quark phase and /0^(r) = p p (r) — p e {r) in the hadron phase. The other 
equations of motion give only the expressions for the chemical potentials, 

^kin+str 



Pi 



NiV(r), 



(54) 



where E kin+ 

str 

expressed as 



E 



d 3 rei + Ej$. Then the quark chemical potentials are 



Pv 



2a r \ 2 i 2 
l + ^Wr)--nr) 



2a r \ 2 i 1 

^ = (i + -^J^iw + ^w 



M . 2a c 



PF S (r) 



eFs(r) 



In 



e Fs (r) +PF S (r) 



(55) 
(56) 

). 
(57) 
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with £Fs(r) = \ffn1 + p| s (r), whereas the chemical potentials of the nucleons 
and electrons are 

_ PFn 

~~ h Cbind 



2m po 
6 



/ Pn(r) + p p (r) _\ Kq ( 
\ Po J 9 V Po 



, ?r p»(r)+p p (r) 

T ^sat '-'sat 1,0° J 

PO 

fx p = Mn „%(£) + g|^)_45 (p B -2, p (r)) 2 _ 
2m 2m po 

p e = (3^ 2 p e (r))^+^(r)- (60) 

We solve these equations of motion under the GC. The important point here is 
that the Coulomb potential V(r) is included properly in each expression. The 
Coulomb potential is a function of the charged-particle densities, and, in turn, 
these densities are functions of the Coulomb potential. As a result, the Poisson 
equation becomes highly non-linear. 

Note that we must now determine eight variables, i.e., six chemical poten- 
tials, fx u , Pd, Pa, Pp, Pn, Pe, and the radii R and i?w- First, we fix R and i?w- 
Here we have four conditions due to the beta equilibrium, expressed by (|52|l . 
Therefore, once the two chemical potentials pb and fi e are given, we can deter- 
mine the other four chemical potentials, p u , pd, p s and p p . Next, we determine 
p e by the global charge neutrality condition: 

d 3 rpl + [ d 3 rpf h = 0. (61) 



ch 

V Q JVn 

The pressure coming from the surface tension 5 is given by 

Then, we find the optimal value of R (i?w is fixed, and therefore fy changes 
with R) by using one of the GC, 

P Q =P H +P T . (63) 

The pressure in each phase, P^ or P H , is given by the thermodynamic relation 
pQ(H) _ _ Qq(h)/Vq(h), where r2q( H ) is the thermodynamic potential in each 
phase and given by adding an electron and the Coulomb interaction contribu- 
tions to ^q(h) m Eqs. I|47|) and (|49|l . Finally, we determine i?w by minimizing 

5 Basically the pressure is locally defined in terms of the Gauss principal curvatures, n\ and 
K2, on the surface: P T = r (ki + Ata) • For the special shapes of sphere, cylinder and plane, Ri 
is simply constant to give (k\, k^) = 1/-R), (V^i 0) an< ^ (0| 0), respectively. 
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the thermodynamic potential. Therefore, once /ib is given, all the other \ii 
(i = u, d, s,p, e), along with R and i?w, can be obtained. 

Note that we keep the GC throughout the numerical procedure. We see 
below how the mixed phase would be changed by including the finite-size effects, 
keeping the GC completely. Although the MC is not strictly correct, our results 
exhibit behavior similar to that in the case of the MC, as a result of including 
the finite-size effects. 

In the numerical calculation, every point inside a cell is represented by a 
grid point (the number of grid points is A gri d w 100). The equations of motion 
are solved by using a relaxation method for a given baryon-number chemical 
potential under the constraint of global charge neutrality. 



7 Charge screening in the hadron-quark mixed 
phase 

7.1 Bulk calculation 

First of all we depict the deconfinement transition between uniform hadron 
matter and quark matter in Fig. by using Eqs. (|47|) . 149fl and (|51|l . Then we 
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Figure 21: Energy density e for uniform hadron matter and quark matter. Uni- 
form quark matter is energetically favorable in high-density region, while hadron 
matter in low-density region. 



can see that it exhibits the first-order phase transition. To describe the mixed 
phase a bulk calculation proceeds like this. Considering two semi- infinite matter 
separated by a sharp boundary and introducing the volume fraction of quark 
matter /y, we apply the global charge neutrality condition, 

fvpl + (1 " fv) P l = 0. (64) 

Assuming the chemical equilibrium among the particles, there are left two in- 
dependent quantities to be determined, the electron chemical potential fi e and 
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the baryon- number chemical potential fiB = fi>n- Finally imposing the pressure 
balance condition between two phases, 



Ph = Pa 



(65) 



we can determine fi e and fy for a given ixb or baryon- number density. We show 
the EOS in this calculation in Fig. (|22|) . Compared with the EOS given by the 
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Figure 22: Pressure of uniform matter, that given by the bulk calculation with 
GC (Bulk Gibbs) and that given by MC (Maxwell). 



Maxwell construction, where the local charge neutrality is implicitly assumed 
instead of Eq. (|B"£}l , we can see that the density regime of the mixed phase is 
considerably widened and pressure is no more constant there. 

The above picture of the mixed phase is too simple because it should have 
structures by the surface tension and the Coulomb interaction, which are called 
the finite-size-effects. Actually Heiselberg et al. demonstrated that the finite- 
size effects may energetically disfavor the mixed phase and thereby its density 
regime is largely limited In their calculation, however, the non-uniform 

structure of the mixed phase is always stable due to an inconsistent inclusion 
of the Coulomb interaction. We shall see that the rearrangement of the charge 
density and the charge screening of the Coulomb interaction together may cause 
a mechanical instability of the structured mixed phase. As another defect they 
use the pressure balance condition before introducing the finite-size effects; they 
use the value of fy given by the above simple calculation. 

7.2 Mechanical instability 

Some analytic discussions have been given in Refs. [6111621 only by taking the 
leading contribution of the Coulomb potential in the Poisson equation. We shall 
see in the next section the numerical results about the properties of the mixed 
phase and EOS, but it should be useful to figure out some qualitative aspects 
by analytic discussions. In fact we can clearly see a mechanical instability of 
the structured mixed phase due to the charge screening effect. 
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If we expand the charge density in 6V(r) = V(r) — V Te { around the reference 
value V IC [, which is also gauge dependent, the Poisson equation (|53|l renders 

A6V a (r) = 4Tre 2 p ch - a (V a (r) = V&) + {n a {V a {r) = V£ f )f8V a (r) + ■■■, (66) 

in each phase a = Q or H with the Debye screening parameter, 



(« a (V«(r) = KSf)) 2 = 47re 



dV 



4 nQ?Q« 



y°( r )=y° 

(67) 

where we used the gauge-invariance relation, 

Then we calculate contribution to the thermodynamic potential (effective en- 
ergy) of the cell up to 0(5V a (r)) 2 . The Coulomb interaction energy can be 
written by way of the Poisson equation l|66|l as 

E v = I / / / ^^*r. («») 



j Vq JVn Jv Q 87re 2 J Va 8ne 2 

Besides the terms given by ()69(l . there are another contributions arising from ef- 
fects associated with the inhomogeneity of the electric potential profile, through 
implicit dependence of the particle densities on V^ ,H (r). We will call them 
"correlation terms", uj" m = e|J in+str — nfpf. Taking e£ in+str and pf as func- 
tions of V a (r) we expand them around the reference value V rc { to obtain the 
corresponding correlation contribution to the thermodynamic potential f2 cor = 
Jvq d 3 r^ OI + J Vh d 3 ru? m . 



cor — e kin+str(Pi*(Kef)) MiVfOref) P^'" (Kef )Kef 

V r g f AV-"(r) | (^(V r a e{ )) 2 (5V a (r)) 2 | 
Aire 2 87re 2 



(70) 



where we also used Eqs. (|Blj|) and l|67|). In general V® { ^ and they may 
depend on the droplet size. Their proper choice should provide appropriate 
convergence of the above expansion in SV(r). Taking 

K?f = Kef = Kef = const (71) 

we find 

(««(v; ef )) 2 (y«(r)-y ref ) 2 



87re 2 



+ const, (72) 



and one may count the potential from the corresponding constant value. Here 
we also took into account that the term V Y " { AV a (r) / (Aire 2 ) does not contribute 
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to fl according to the boundary conditions in the droplet center, at the droplet 
boundary (at zero surface charge) , and at the boundary of the Wigner-Seitz cell. 

The Poisson equation (jfjj))l in the droplet phase can be easily solved. 6 For 
r < R, the charge density p^ h can be written as 



2/4(FQ(r)- Me ) A , „ ( (VQ(r) - fi e )\\ , /W 1 



1 + 



9-7T 2 ^ y [Lb ) ) 67T 2 

by way of Eq. (|55|l to find 



(73) 



V Q (r) - a e = sinh( K Q r) + C/ n Q , (74) 

with an arbitrary constant V^. For the Debye parameter re^ and for the con- 
stant we obtain: 



<« u > J -=ST-rJ- W = (75) 

Thus, the value is rather small and the main contribution to V^(r) comes 
from the first term in (|74|l . Note that solution l|74[) is independent of the refer- 
ence value V^f in this case, cf. JB7JI, since p® h in (T7S|) is the linear function of 
VQ(r) — in the approximation used. 

For r > R, expanding the charge density p^(r) around a reference value, 

pi (r) ~ Pp (y H (r) = O + 5p p (r) - p e (V H (r) = Vg) - Sp e (r), (76) 

we find up to linear order 

Sp p (r) ~ Co-H^W - KS), <fp e (r) = ^- ^cf) 2 (t/H(r) _ yH f))(77) 

7T 

C ° = W' ^ = ( 37r %(^ H ( r )=^ef)) 1/3 , 

where ^22 1 an d | ^4 | are the corresponding matrix element and the determinant 
of the matrix, 

K ~ s Pi s Pj ■ (78) 

Then the Poisson equation with the boundary condition 

VV H (v)\ r=Rw = (79) 

6 We only consider the quark-droplet phase here for simplicity. See Refs. 61 62 90 for 
other cases. 
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yields 



F H (r) - Me = F H ^ cosh ( K H (r - i? w )) (1 - 5) + C/ H , (80) 



8 = tanh (K H (i? w - r)) /(n H R w ), 
with an arbitrary constant Vg H , where the constant is given by 



^0 + Me T^hT^ 1" ^ref • l 8i J 



The charge screening in the external region is determined by the Debye 
parameter 

where the second term is the contribution of the proton screening. Taking 
p§ = 1.5po, Me,Gibbs — 170 MeV, fig = /it n ~ 1020 MeV, a c ~ 0.4, we estimate 
typical Debye screening lengths as A§ = 1/k^ ~ 3.4/m w , and A§ = 1/k h ~ 
4.2/m ff , whereas one would have ~ 8.5/tott, if the proton contribution to 
the screening (|S2*|) was absent (C^ 1 = 0). With the estimate A§ ~ we 
get that K H i?w > 1 f° r the droplets with the radii R > (2/y) 1 / 3 • 3.3/m w . 

In Fig. E3] we demonstrate the radius (i?) dependence of the total ther- 
modynamic potential per droplet volume for the case of spherical droplets, 
<5u3tot/A) = (e^ + £y + ^Sr + ^?or + es)/A), given by the sum of partial con- 
tributions, where tilde denotes each quantity scaled by the droplet volume 
V = AttR 3 /3 and /3o is a typical quantity with the dimension of the energy 
density H3 EH El, 



3(l/ H -l/ Q ) 2 (^) 2 
8?re 2 



0o = V " V ; ■ (83) 



Preparing some wide range for the value of the surface tension parameter r or 
ft = 3k q t//3 ED El EDI, we present two cases of f v , f v = 0.01,0.5. The 
dotted green line shows the asymptotic value of (5a3 tot //3 in the limit £ — > oo. 

The label "C" is given for reference to show the previous non self-consistent 
case, where the Coulomb potential is not screened (see Fig. EJ). We can see 
that only in the limit of fv <IC 1 and R <C A§, we are able to recover this 
case. The "e.m." curve shows the partial contribution to the thermodynamic 
potential, ec.m./ft) = (eV + es)/ft), ignoring correlation terms. Comparing these 
curves we can see how the screening effect changes the thermodynamic potential: 
we can see that the minima at the "e.m." curves disappear already at ft > 
0.03, corresponding to unphysically small r ~ several MeV/fm 2 . However, the 
correlation energy gives a sizable contribution to allow the minimum for larger 
value of t. Consequently, the minimum totally disappears between ft =0.1 and 
ft = 0.5, which may be interpreted as 10 < r < 50 MeV/fm 2 in this calculation. 
Thus we have seen a mechanical instability of the droplet for the medium values 
of r, which might be in the physically meaningful range. 
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Figure 23: Dimensionless value of the thermodynamic potential per droplet 
volume. Solid lines are given for fy = 0.5 and dashed lines for fy = 1/100. 
The ratio of the screening lengths of two phases, otQ — A§/A^ , is fixed as one. 
£ is a dimensionless radius of the droplet, £ = i?/A§, with A§ ~ 5 fm in this 
calculation. See text for further details. 



7.3 Charge screening in quark pasta structures 

We display the thermodynamic potential in Figs. l24l and l25l In uniform matter, 
the hadron phase is thermodynamically favorable for /ib < 1225 MeV, and the 
quark phase for /^b > 1225 MeV. Therefore we plot Su>, the difference between 
the thermodynamic potential densities of the mixed phase and each phase of 
uniform matter: 

, / "total - ^ H niform MB < 1225 MeV, 

\ "total - "£ niform MB > 1225 MeV. {m > 

Here, "total = ^totai/Vw, etc. There we also depict two results for compari- 
son: one is obtained from the bulk calculation, where the finite-size effects are 
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Figure 24: Difference between the 
thermodynamic potential densities as 
a function of baryon-number chemical 
potential for r — 40 MeV/fm 2 . 
If bto is negative, the mixed phase 
is the thermodynamically favorable 
state. The MC determines one point 
of the phase transition in uniform mat- 
ter, denoted as a circle in the figure. 



Figure 25: The same as Fig. [23] for 
r — 60 MeV/fm 2 . The region of nega- 
tive 5u> is narrower than in the r = 40 
MeV/fm 2 case. 



completely discarded (cf. Fig. I22|l. and the other is the thermodynamic poten- 
tial obtained using a perturbative treatment of the Coulomb interaction, which 
is denoted by "no Coulomb" . We have employed the similar procedure as in 
Refs. [SS]E|[3HI. In this procedure, by discarding the Coulomb potential V(r) 
to solve the equations of motion Eqs. J5SJ|-||B2J, each density is determined to 
be constant in each phase for given R and Ryj. Using the chemical equilib- 
rium relation 1|52|) and the charge neutrality condition l|t . we can determine 
its value. The Coulomb interaction energy (|42|) can be separately evaluated by 
using the constant densities, and the total thermodynamic potential is obtained 
by adding it. The remaining procedure is the same as that described in Sec. 
6.2: we determine R using the pressure balance relation Eq. I|63|> and the cell 
size (i?w) to minimize the total thermodynamic potential. 7 

We can see the screening effects by comparing these "no Coulomb" results 
with the self-consistent treatment denoted by "screening" . The quantity Su 
derived using the MC appears as a point denoted by a circle in Figs. [2] and 
1251 where only the two conditions = P H and /ig = /ig are satisfied. The 
mixed phase derived with the bulk Gibbs calculation appears over a wide range 
of values of the /ie • Therefore, a narrowing the region of the mixed phase signals 
that the properties of the mixed phase have become closer to those of the MC. 
It is clearly seen that Wtotal becomes close to that given by the MC due to the 
finite-size effects, the effects of the surface tension and the Coulomb interaction. 
Comparing the result of the self-consistent calculation with that in the case of 
"no Coulomb" , we can see that the change in energy caused by the screening 

7 Note that they also discarded the size dependence of the particle densities and the surface 
tension for the pressure balance condition in the bulk calculations 1551 1571 1551 . 
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Figure 26: Density profiles in 
the droplet phase obtained with "no 
Coulomb" for ^ B = 1189 MeV and 
r = 60 MeV/fm 2 . It is seen that 
they are uniform in each phase. Here, 
R = 7.2 fm and R w = 12.8 fin. 



Figure 27: Density profiles and the 
Coulomb potential derived using the 
self-consistent calculation for the same 
parameter set as in Fig. |2E1 Here 
R = 7.7 fm and i? w = 18.9 fm. 



effect is not so large, but still the same order of magnitude as that caused by 
the surface effect. The charge screening effect decreases the Coulomb energy. 
However, it increases the electron energy instead. 

Our parameters are the same as those of Heiselberg et al. |21 , and the mixed 
phase disappears if r > 90 MeV/fm 2 in the case of "no Coulomb". However, 
in the case of screening, it is rather difficult to estimate the precise value of r 
at which the mixed phase disappears, because a mechanical instability would 
appear for values of r that are not too large I^Ul 1611 162j . Nevertheless we 
could infer that the mixed phase would disappear for r = 70 - 80 MeV/fm 2 by 
referring to the difference between the cases of 40 MeV/fm 2 fFig. I24|) and 60 
MeV/fm 2 (Fig.HSJ. 

If the surface tension becomes stronger, the relative importance of the screen- 
ing effect becomes smaller, and the effect of the surface tension becomes more 
dominant, as seen in Figs. 1241 and 1251 Although the charge screening does not 
have so large effect on the bulk properties of matter, we shall see that its effect 
is significant for the charged particles in the mixed phase, and this brings about 
a significant change of the properties of the mixed phase. 

The charge screening effect induces a rearrangement of the charged parti- 
cles. 8 We can see this screening effect by comparing Fig. with Fig. [Z7\ The 
quark phase is negatively charged, and the hadron phase is positively charged. 
The negatively charged particles in the quark phase, such as d, s and e, and 
the positively charged particle in the hadron phase, p, are attracted toward the 
boundary. In contrast, the positively charged particle in the quark phase, if, 
and the negatively charged particle in the hadron phase, e, are repelled from 
the boundary. 

8 In this context, it would be interesting to refer the paper by Heiselberg who studied 
the screening effect on a quark droplet (strangelet) in the vacuum, and also suggested the 
importance of the rearrangement of the charged particle densities. 
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Figure 28: Local charge densities in the case of "no Coulomb" and the case of 
the self-consistent calculation with the screening effect. For "no Coulomb" , the 
charge density of each phase is constant over the region. The absolute value of 
the charge density is larger than that derived from the self-consistent calculation 
in each phase. In the hadron phase, the charge density almost vanishes near the 
cell boundary, r = Ryj. 



The charge screening effect also reduces the net charge in each phase. In 
Fig. 1281 we display the local charge densities of the two cases considered in 
Figs. EE1 and |23 The change of the number of charged particles due to the 
screening is as follows: in the quark phase, the numbers of d and s quarks 
and electrons decrease, while the number of u quarks increases. In the hadron 
phase, contrastingly, the number of protons should decrease and the number 
of electrons should increase. Consequently, the local charge decreases in both 
phases. In Fig. EHl we can see that the core region of the droplet tends to 
be charge neutral, and near the boundary of the Wigner-Seitz cell it is almost 
charge neutral. 

In Figs. and |20l we present radii of the lump and cell as functions of 
density. The R dependence of the total thermodynamic potential comes from 
the contributions of the surface tension and the Coulomb interaction. The op- 
timal radius representing the minimum of the thermodynamic potential is then 
determined by the balance between these two contributions. If the Coulomb 
interaction energy is suppressed, the minimum of the thermodynamic potential 
is shifted in the direction of larger radius. As a result, the size of the embedded 
phase R and the cell size i?w become large. In the previous section we have seen 
that the minimum disappears for a large value of the surface tension parameter, 
and hence the structure becomes mechanically unstable. We cannot show this 
directly in our framework, because such unstable solutions are automatically 
excluded during the numerical procedure, although we can see this tendency in 
Figs. 1291 and 1301 R and i?w become larger through the screening effect. 

We also see the relation between the size of the geometrical structure and the 
Debye screening length. The Debye screening length appears in the linearized 
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Figure 29: Lump and cell radii de- 
rived in the "no Coulomb" calculation. 
The Debye screening length is also de- 
picted for comparison. R is repre- 
sented by the thick, solid line and i?w 
by the thick, dashed line. We can see 
that the size of the structure becomes 
less than the Debye screening length. 



Figure 30: The same quantities as 
in Fig. |2H1 derived using the self- 
consistent calculation with the screen- 
ing effect. The size of the structure 
here is larger than that obtained with 
"no Coulomb" , and exceeds the Debye 
screening length. 





Figure 31: Schematic graphs of the droplet size and the Debye screening length. 
The right figure corresponds to the case of the self-consistent calculation with 
the screening effect, and left figure corresponds to the case of "no Coulomb". 
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Poisson equation and is given by 



™-^x>(^), (^^.(^). 

(85) 

where (p^ h ) is the average density in the quark phase, (pp h ) is the average proton 
number density in the hadron phase, and (pf 1 ) is the average electron number 
density inside the cell. The Debye screening length gives a rough measure of 
the screening effect: at a distance larger than the Debye screening length, the 
Coulomb interaction is effectively suppressed. 

In Fig. we present sizes of the geometrical structure in the case of "no 
Coulomb". If we ignore the screening effect, the size of the embedded phase is 
comparable to or smaller than the corresponding quark Debye screening length, 
Ap, for a droplet, rod and slab, or the proton Debye screening length, Aq, for 
a tube and bubble (Fig. EHI) . This may imply that the Debye screening is not 
so important. Actually, many authors have ignored the screening effect for the 
reason elucidated by this argument P^HSHl- In Fig. [301 however, we see that 
the size of the embedded phase can be larger than A^ (Fig. I^T|) in the self- 
consistent calculation. We can also see a similar situation concerning i?w and 
Ap. This means that the screening has important effects in this mixed phase. 
We cannot expect such an effect without solving the Poisson equation, because 
of the non-linearity. 9 

We plot the EOS in Figs. 1521 and!531 It is seen that the pressure of the mixed 
phase becomes similar to that found using the MC. Because a region of local 
charge neutrality appears due to the screening effect (see Fig. EH}, its properties 
correspond to those derived using the MC. 

8 Summary and concluding remarks 

We have investigated three forms of non-uniform hadron matter at different 
density ranges, which correspond to the mixed phases in the first-order phase 
transitions. As has been discussed by many authors, such non-uniform struc- 
tures appear as "pasta" structures, i.e. droplet, rod, slab, tube and bubble. By 
means of mean-field approach under the Wigner-Seitz cell approximation, we 
have calculated the ground-state density profiles of matter. The size and the 
dimensionality of the cell were chosen to their optimal values for a given baryon 
density. Thus we could determine optimal structure of matter along the change 
of density. 

First we have studied the pasta structures in low-density nucleon matter. 
In this study we have used the relativistic mean-fields (RMFs) as an effective 

9 Within the linear approximation for the charge screening, Iida ans Sato have studied the 
screening effect on quark droplets in a different context [92) , In their case, the electron number 
density is much larger than in our case, because they considered two-flavor quark matter. 
Accordingly, the electron screening effect may be dominant, in contrast to the situation in our 
case. 
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Figure 32: Pressure as a function of 
the baryon-number density obtained 
from the "no Coulomb" calculation for 
t = 40 MeV/fm 2 . The results found 
with the bulk Gibbs calculation and 
the MC are also presented for compar- 
ison. 




Figure 33: The same quantities as 
in Fig. computed using the self- 
consistent calculation with the screen- 
ing effect. 



interaction. The parameter set we have used gives the incompressibility of 
matter K — 240 MeV and the symmetry energy coefficient of 32.5 MeV at 
the saturation density of symmetric nuclear matter, and can well reproduce 
the bulk properties of finite systems (nuclei). In nucleon matter with a fixed 
proton mixing ratio, which corresponds to supernova matter in the collapsing 
stage, we have observed the nuclear pasta: with increase of density the favorable 
structure changes from nuclear droplet, rod, slab, tube, and to bubble. Near the 
saturation density (normal nuclear density) , nuclear pasta dissolved to uniform 
matter. By the appearance of such non-uniform structure, the equation of state 
of the system becomes significantly soft (the energy gets lower up to 15 MeV/ A). 
In the case of the beta equilibrium (neutron star matter), we did not observe 
the full pasta structures but only proton droplets appeared in the neutron sea. 

Our second target was the kaon condensation at around three times the 
normal nuclear density, which is relevant to the inner core of neutron stars. 
Using the same RMF model, we have added the kaonic degree of freedom and 
could naturally describe both normal nucleon matter and kaonic matter. The 
surface tension between nucleon matter and kaonic matter is automatically in- 
cluded in the framework. With this model we observed the appearance of pasta 
structures: kaonic droplet, kaonic rod, etc. 

An important advantage of our framework is the consistent treatment of the 
Coulomb interaction and charged particle distribution: we have included the 
Coulomb potential in the thermodynamic potential in a gauge-independent way 
and solved the Poisson equation consistently with the charged particle densities. 
If we take uniform charge densities as in the bulk calculations, we overestimate 
the potential energy due to the neglect of the charge screening effect. In fact, 
we have seen the importance of the charge screening for the properties of the 
mixed phases; the region of the mixed phase is largely modified, and it may 
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cause a mechanical instability of the geometrical structure in some cases. As an 
important consequence of the charge screening effect, the EOS becomes more 
or less similar to that given by the Maxwell construction which is irrelevant for 
multiple particle composition. Actually the electron chemical potential should 
be different in the absence of the Coulomb potential because of the different elec- 
tron number in two phases. However, we have emphasized that charge chemical 
potential is only defined after the gauge fixing of the Coulomb potential, so 
that number densities are not solely given by the chemical potential but by the 
linear combination of the chemical potential and the Coulomb potential. The 
different strength of the Coulomb potential in two phases can give the different 
particle densities, which can resolve an apparent contradiction in the Maxwell 
construction. 

So far, we have explored only nuclear matter in the ground state (zero tem- 
perature). However, supernova matter, for example, should be as hot as several 
MeV. ft is naturally expected that the pasta structure will be modified and dis- 
solved at around the critical temperature [7J . Our model is also applicable for 
not-very-high temperature by modifying the baryon momentum distribution. 
The structural phase diagram in density-temperature plane should be useful for 
the discussion of supernova explosion. 

We have used the meson exchange model (MEM) to describe the kaon- 
nucleon interaction. If we use the chiral model instead, it gives a nonlinear 
Lagrangian including the kaon self-interaction. The thermodynamic inconsis- 
tency of the chiral model has been shown within the mean-field approximation 
|2fill^2llM| : a bulk calculation of the mixed phase consisting of two semi-infinite 
matter have concluded that we can not construct the mixed phase satisfying the 
Gibbs conditions. It may be caused by the non-linearity of the kaon field. So, 
it would be worth studying the chiral model with our framework to see whether 
the chiral model is really ill defined in a thermodynamic sense. 

We have seen that our framework based on RMF can well reproduce the bulk 
properties of finite nuclei as well as nuclear matter. It should be interesting to 
describe kaonic nuclei [941 1951 1551 I97| as a direct application of our framework. 
Then we can get a consistent description of kaonic nuclei by taking into account 
the Coulomb interaction as well as the K-N interaction 98 . 

Then our third subject was the mixed phase of hadron-quark phase transi- 
tion. We have used a rather simple model for this study. Nowadays there have 
been many studies about the deconfincment transition and its implications on 
compact stars 0E1G>5 , 56 , 57 and relativistic heavy- ion collisions |47U52ll5lfll54| . 
Let us consider some implications of these results for the structure of hybrid 
stars. Glendenning and Pei considered some geometrical structures in the 
mixed phase and conjectured that "crystalline structures" (quark pastas) would 
appear in the core region of a hybrid star by using the bulk Gibbs calculation. 
His results predict that the mixed phase should appear for several kilometers. 
However, we can say that the region of SMP should be narrow in the /ie space, 
and the EOS is more similar to that obtained with the MC, due to the finite-size 
effects. 

We have considered the deconfinement transition from non-strange nuclear 
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matter for simplicity, while realistically we must take into account hyperons. 
Many people believe that hyperons will mix in nuclear matter at several times 
of the saturation density |15] EM HUD EH EH EH EE3- In such a case 
we must consider the mixed phase by taking into account hyperons as well. 
One of the interesting consequences may be related to the properties of hybrid 
stars: actually it has been shown that the admixture of hyperons considerably 
softens EOS, while the deconfinement transition stiffens EOS again at high- 
density region |1U4I I105| . If the stiffening works well, it becomes rather easy to 
construct hybrid stars. For quark matter we should consider the possibilities of 
color superconductivity an d ferromagnetism HUB] . In particular, 

magnetic properties of quark matter might be related to the strong magnetic 
field observed in compact stars [1071 11081 1109| . 

Our framework may be applicable for other subjects. One of the interesting 
subjects may be color superconductivity. There have been many works about 
color superconductivity, where various first-order phase transitions have been 
suggested between different types of superconducting phases |11(JI lllll I112| . 
Color neutrality as well as charge neutrality should be considered in that case. 
Accordingly 3rd and 8th gluons should play an important role for global color 
neutrality. 



Acknowledgments 

The authors thank M. Alford, M. Baldo, D. Blaschke, A. Bonasera, K. Iida, 
A. Iwamoto, N. Iwamoto, E. Kolomeitsev, S. Kubis, U. Lombardo, Tomoyuki 
Maruyama, T. Muto, S. Nishizaki, K. Oyamatsu, C. J. Pethick, S. Reddy, H.- 
J. Schulze, T. Takatsuka, R. Tamagaki, T. Tanigawa, S. Tsuruta, D. N. Voskre- 
sensky, G. Watanabe, F. Weber and M. Yasuhira for their interests and discus- 
sions. This work is partially supported by the Grant-in- Aid for the 21st Century 
COE "Center for the Diversity and Universality in Physics" from the Ministry 
of Education, Culture, Sports, Science and Technology of Japan. It is also par- 
tially supported by the Japanese Grant-in- Aid for Scientific Research Fund of 
the Ministry of Education, Culture, Sports, Science and Technology (13640282, 
16540246). 



APPENDIX: Kaon interactions within the chiral 
model 



We have used the meson exchange model (MEM) to describe kaon condensation. 
When we use the chiral model instead, the thermodynamic potential is given as 
follows: 



n = n, 
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where the kaon field is defined as 
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Then the equations of motion can be easily written down, 
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We can see how the thermodynamic potential i|86[) or the equations of motion 
Ij89(l can recover the previous formulae used in the studies of kaon condensation 
in uniform matter 26 , 32 j. The K-N interaction terms are easily extracted from 
Eq. ||SSJ|. Discarding V, the kaon source terms for the mean-fields a, uj,p and 
the nonlinear potential for er, U(a), we have 
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and 

v gioNg^K , s g P Ngu>K , \ 

X = 5 (p n +p p )-\ - (Pp-pn), 

mi m 2 p 

*2 2 o g<jNg<jK I s . s \ rm\ 

m*K = m z K -2m K — {p n +Pp), (91) 

for soft kaons. If we impose the following relations among the coupling con- 
stants: 
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X ° = 772^» + 2 ^p)' 



(92) 

we can recover the K-N interaction terms dictated by chiral symmetry, 

1 

mf = m 2 K -^f(p s n + P s p ). (93) 

Thus we can derive the previous formulae for kaon condensation within the 
chiral model |261 15*2^ 

To get MEM, we linearize Eqs. (|56*|) and lj%9")l with respect to 9 and further 



add the non-linear terms, X%9 2 and a 2 9 2 , in fi^ . 
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